A1 = readmatrix('NiI2_MT_HEAT.dat')
A2 = readmatrix('NiI2_09T_cool.dat')

x1 = A1(:,3)
y1 = A1(:,59)*130

x2 = A2(:,3)
y2 = A2(:,5)*11

hold on
%scatter(x1,y1,5,[0 0 0.8],'filled')
plot(x1,y1,Color=[0.8 0 0])
plot(x2,y2,Color=[0 0 0.8])
xlabel('Temperature [K]')
ylabel('chi [emu/mol]')
ylim([5.5e-3,8.2e-3])
xlim([10,100])
yticks([6e-3:1e-3:8e-3])
xticks([20:20:100])

set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 200, 160])